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ABSTRACT. Aside from primordial gravitational instability of the cosmological fluid, various mechanisms have 
been proposed to generate large-scale structure at relatively late times, including, e.g., "late-time" cosmological 
phase transitions. In these scenarios, it is envisioned that the universe is nearly homogeneous at the time of 
last scattering and that perturbations grow rapidly sometime after the primordial plasma recombines. On this 
basis, it was suggested that large inhomogeneities could be generated while leaving relatively little imprint on the 
cosmic microwave background (MBR) anisotropy. In this paper, we calculate the minimal anisotropics possible 
in any "late-time" scenario for structure formation, given the level of inhomogeneity observed at present. Since 
the growth of the inhomogeneity involves time-varying gravitational fields, these scenarios inevitably generate 
significant MBR anisotropy via the Sachs- Wolfe effect. Moreover, we show that the large-angle MBR anisotropy 
produced by the rapid post-recombination growth of inhomogeneity is generally greater than that produced by the 
same inhomogeneity grown via gravitational instability. In "realistic" scenarios one can decrease the anisotropy 
compared to models with primordial adiabatic fluctuations, but only on very small angular scales. The value 
of any particular measure of the anisotropy can be made small in late-time models, but only by making the 
time-dependence of the gravitational field sufficiently "pathological" . 



I. Introduction 

Soon after the discovery of the Microwave Background Radiation (MBR), it was noted that mea- 
surements of differences in the MBR temperature in different directions (anisotropy) would provide 
a sensitive probe of large-scale density inhomogeneities in the universe (Sachs and Wolfe 1967, Rees 
and Sciama 1968). Recently, the COBE satellite has discovered MBR anisotropy on large angular 
scales (> 10°) (Smoot, et al. 1992), and experiments on smaller angular scales have seen signals 
which may also turn out to be anisotropy in the MBR (Gaier, et al. 1992, Devlin, et al. 1992, Meyer, 
et al. 1991). While it is clear that large scale density perturbations in the universe will induce MBR 
anisotropies, the relationship between the anisotropics and the perturbations depends on how the in- 
homogeneities are produced. The usual assumption is that the density perturbations are primordial, 
i.e., produced long before recombination, and evolve gravitationally in a relatively simple cosmologi- 
cal fluid composed of, e.g., photons, neutrinos, baryons, and dark matter. If the cosmological matter 



does have such a simple equation of state, then the inhomogeneities must be primordial, since they 
will not arise spontaneously in such a fluid. The relation between density inhomogeneities and tem- 
perature anisotropies for such primordial perturbations is fairly well understood (Sachs and Wolfe 
1967; for a recent introduction, see Peebles 1993), and the implications of the COBE detection for 
primordial perturbations have been extensively studied (e.g., Wright et al. 1992, Efstathiou, Bond, 
and White 1992). In the simplest case, a spatially flat, Einstein-de Sitter universe with = 1, the 
MBR anisotropy gives an imprint of conditions at recombination, when the MBR last scattered; we 
will call these "primordial" anisotropies. 

The other possibility to consider is that the dynamics of some component of the matter in the 
universe is much more complex than that of a simple fluid, and is able to induce perturbations even at 
very late times. One class of such models is that of topological defects, such as cosmic strings, textures, 
or global monopoles. The dynamics of the defects is nontrivial, and as they move around, they 
induce perturbations in the other matter components (baryons, photons, etc.) via their gravitational 
attraction. In these models, perturbations are produced both before and after recombination. The 
induced MBR temperature fluctuations are a mixture of the classical "primordial" anisotropies and 
anisotropies produced after recombination (Stebbins 1988, Turok and Spergel 1990, Bouchet, Bennett, 
and Stebbins 1988, Bennett, Stebbins, and Bouchet 1992, Bennett and Rhie 1992). 

Another class of models are ones in which perturbations are generated primarily after recom- 
bination, or more specifically, in which the inhomogeneity in the gravitational potential increases 
significantly after recombination. These scenarios include "late-time" phase transitions (Wasserman 
1986, Hill, Schramm, and Fry 1989, Hill, Schramm, and Widrow 1991, Press, Ryden, and Spergel 
1990, Fuller and Schramm 1992, Frieman, Hill, and Watkins 1992) involving non-trivial scalar field 
dynamics. In this scenario, one could start with an essentially homogeneous universe at the epoch 
of last scattering and thus avoid all "primordial" anisotropy. Consequently, it was thought, such 
models could generate the observed large-scale structure with very small imprint on the MBR, and 
late phase transitions were posed as alternatives to the standard primordial gravitational instability 
scenario, which, even after COBE, appears to be on the edge of producing an excessive small-scale 
MBR anisotropy. However, one does not completely avoid MBR anisotropies in the late-time scenario: 
in this case, both density perturbations and MBR anisotropies are produced after recombination. If 
the universe is homogeneous at last scattering, then gravitational field perturbations must subse- 
quently grow from zero to their present value in order to account for the observed structure. This 
time- varying gravitational field will induce MBR anisotropies which may not be very small compared 
with the "primordial" anisotropies produced in the primordial instability scenario. In other words, 
MBR anisotropies are an inevitable consequence of the existence of density perturbations today. 

In this paper, we calculate the minimal MBR anisotropies associated with scenarios for late 
structure formation such as a late-time phase transition. More specifically, we determine the minimal 
anisotropies implied by the boundary conditions of zero inhomogeneity at recombination and a fixed 
present amplitude of density inhomogeneity inferred from redshift surveys. We do not deal with any 
particular model of such a phase transition, so we give no predictions for the anisotropies that might 
be expected in a completely 'realistic' scenario. Rather we set a firm lower limit on how small these 
anisotropies can be, given the observed level of large-scale structure. As we shall see, this lower limit 
is not particularly small when compared to what is expected in scenarios with primordial density 
inhomogeneities. Moreover, when reasonable smoothness conditions are placed on the evolution of 
the gravitational potential, the large-angle anisotropy in the late-time scenario is generally larger than 
that for primordial adiabatic fluctuations with the same present level of clustering. 

Several estimates of the induced anisotropy have previously been calculated for specific late-time 
scenarios (Zel'dovich, Kobzarev, and Okun 1974, Stebbins and Turner 1989, Turner, Watkins, and 
Widrow 1991). These estimates suggest that, in models involving domain walls, the anisotropy may 
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be significantly larger than in the primordial gravitational instability picture. However, methods to 
fix up the domain wall problem have been suggested (Massarotti 1991, Massarotti and Quashnock 
1992), and late-time transitions without domain walls have also been investigated (Press, Ryden, and 
Spergel 1990, Frieman, Hill, and Watkins 1992). This flux of theoretical developments has motivated 
the model-independent approach we adopt here. 

In §11, we set up the problem of minimizing the temperature correlation function C(a). We 
analytically study the minimization of the rms variance in the temperature fluctuation, C(0), first 
assuming the time-dependence of the gravitational potential is continuous and bounded, and then 
allowing it to become unbounded. In the latter case, the anisotropy is reduced. This illustrates our 
general result: if the gravitational potential is constrained to be a smoothly varying function of time, 
the anisotropy is substantially larger than if the potential is unconstrained. In particular, the large- 
angle anisotropy for a sufficiently smoothly varying potential is generally larger than for primordial 
adiabatic perturbations, while the anisotropy for a pathological potential function can be smaller. 
We also point out the differences expected in the angular dependence of the anisotropy between the 
primordial and late-time scenarios. In §111, we discuss minimization of the anisotropy for experiments 
with finite beam-width, and solve the problem numerically for several beam configurations. In §IV, we 
use analytic fits to the power spectrum of large-scale clustering suggested by recent redshift surveys 
to normalize the results and make estimates of the minimal anisotropy for different experimental 
beamwidths. We conclude in §V. The details of some of the numerical computations are relegated to 
the Appendices. 

(A note on nomenclature: subtracting off the monopole and dipole anisotropy, MBR anisotropies 
on large angular scales can generally be decomposed into two components, one of which can loosely 
be thought of as arising from the gravitational potential at the surface of last scattering, and the 
other as due to the time-dependent gravitational potential along the path of the MBR photons since 
recombination. We are calling these two terms "primordial" and "post-recombination" respectively. 
In the literature, these are sometimes called the "Sachs- Wolfe" and "Rees-Sciama" effects, but they 
are both contained in the Sachs- Wolfe (1967) formalism. Nevertheless, we will sometimes partially 
lapse into this usage as well, and use the terms "primordial" and "usual Sachs- Wolfe" interchange- 
ably. The post-recombination effects we are talking about are rather different than the "Rees-Sciama 
effect": Rees and Sciama (1968) studied the anisotropies produced by spherical growing mode density 
perturbations after recombination, and found differences from the classical Sachs- Wolfe effect only 
due to the non-linearities in the gravitational instability of the matter. The anisotropies we consider 
are linear in the amplitude of the present density inhomogeneities.) 

II. Minimal Models: Formalism and Analytic Results 

As shown by Sachs and Wolfe (1967), a non-uniform gravitational field causes anisotropies 
in the MBR by making differential changes in the energy of photons. In a spatially flat (k = 0) 
Friedmann-Robertson- Walker (FRW) cosmology, which for simplicity we shall assume throughout, 
the perturbed metric can be written = a 2 {r]){ri^ v + h^), where a(rj) is the FRW scale factor, 
Vnv = diag[— 1, 1,1,1] is the Minkowski metric, h^ v is the metric perturbation, and r\ = x° is the 
conformal time, rj = f dt/a. In this geometry, the fractional change in energy of photons moving 
along the null geodesic x^(X) is 

AT 1 f Vi 

-jr = -2j nV (V,o " 2h °^) d V (2-1) 

to first order in the perturbation. In this expression, = dx^/dr/ is the tangent vector of the 
unperturbed photon trajectory, with components n° = 1, n l , where n is a spatial unit vector, r^-nV = 
1, i, j = 1, 2, 3. The derivatives are evaluated at points x^ along the unperturbed photon path. 
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We can decompose a general metric perturbation h fJ/l/ (x fJ ') into scalar, vector, and tensor modes, 
which evolve independently in linear perturbation theory. These modes are orthogonal in the sense 
that the expected value of any quadratic measure of the anisotropy is just a sum of the scalar, 
vector, and tensor components, with no cross terms. Here we only consider quadratic measures of the 
anisotropy. At present, since there is no compelling evidence for the existence of vector and tensor 
metric perturbations in the universe, we can estimate the minimal AT/T consistent with observations 
by considering scalar perturbations alone. Nonzero vector and tensor modes can only increase the 
anisotropy. 

For scalar perturbations, in longitudinal (or conformal Newtonian) gauge, the metric perturba- 
tion takes the form /too = — 2$, hij = —2^>rjij, where $(x, 77) and \I/(x, 77) are gauge-invariant variables. 
(In the notation of Bardeen (1980), $ = $^ and \1/ = — $#.) Substituting into the Sachs- Wolfe in- 
tegral (2.1), setting the metric and temperature fluctuations to zero at recombination, and ignoring 
the boundary term at the observer, which only contributes to the monopole and dipole anisotropy, 
we have 

AT f m 

— = d V (y + $) (2.2) 

where an overdot denotes differentiation with respect to conformal time. This expression is manifestly 
gauge-invariant. For a spatially flat (O = 1) matter-dominated universe, recombination occurs at 
conformal time r\ r = 770/(1 + Zr) 1 / 2 with z r ~ 1100, and the conformal time today is r]Q = 2iT Q ~ 1 = 
6000h _1 Mpc. (We use units in which the speed of light c = 1.) 

If we assume that the universe is dominated by nonrelativistic matter today (or, more generally, 
if the anisotropic stress of the currently dominating matter vanishes, 8T l - ~ <5p, then the present 
boundary condition for the two 'gravitational potentials' $ and \& is $o = ^o- As a result, since the 
two potentials enter the Sachs- Wolfe integral (2.2) identically, for the solution which minimizes AT/T 
they will be equal for all time, $ = (Note that if the universe is currently dominated by matter 
with anisotropic stress, one could violate this assumption. In fact, since only $ contributes to the 
motion of non-relativistic matter, e.g., galaxies, in this case one could conceivably set ^> = for all 
time, reducing the anisotropies we will calculate by a factor of two. However, this would presumably 
require a rather bizarre stress tensor for the dark matter, so we will not consider this possibility 
further.) Putting in the explicit argument of the potential, we thus have 



AT P° 

= 2/ drj^(K -h( Vo -rj),rj), (2.3) 

Jr\ T 



where x Q is the observer position coordinate. In the presence of primordial adiabatic perturbations 
at the surface of last scattering (i.e., what is usually called the Sachs- Wolfe effect), there would also 
be a term (AT/T)gw = 0-/3)®(Vr)- In the case we are considering here, however, there is no metric 
or radiation perturbation initially, and the boundary term at emission (77 = r] r ) is zero. 

For the moment we do not subtract off the contribution to the monopole and dipole contribution 
from (2.3). For experiments that only probe wavenumbers k such that kr\ Q ^> 1, i.e, wavelengths 
much smaller than the present Hubble radius, this subtraction would not make much difference 
anyway. Below, we will consider the rms temperature fluctuation for a realistic experimental beam 
configuration; here, we consider the full temperature autocorrelation function 

/AT AT \ f Vo f v ° 

Clt («) = (— (n) — (m) ) =4/ d V drjf ($(x Q -n (7/0-77), »])$(x -m(»?o-»?'). rf)), 

\ / n-m=cosa Jr\x J r) T 

(2.4) 
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where (• • •) denotes an average over all positions x Q and all directions n, m separated by an angle a. 
Taking the Fourier transform of the potential, 



d 3 k , 

results in 



*(x,v) = I Ts^CMJe**, (2.5) 



(2.6) 



(27T) 

x (e -k.n(, -,) +^ k^m(, -V) )^ , ^il=cosQ( | )(k7??) | ) * (k / 5 ^ )) . 

Performing the appropriate averages: 

(e ix .(k-k') )xo = _L J rf 3 xet (k-k')-x = M! 5 o) (k _ k / } (2 7) 

and (see Appendix A) 



(e- ik - A ^ +ik '^°-^h. A = cosa = h&yfijh - v) 2 + (Vo - v') 2 - 2(Vo - v)(Vo - rf) cosa) (2.8) 
gives 

4 f d 3 k P° p° I 

Clt(c^) = — / — -g / drj dr]'j {kJ(r] o -r]) 2 + (r] - r]') 2 - 2(rj - rj) (r) Q - rf) cos a) 

V J (2ix) 6 J Vi J Vi V ^ 2 _ 9 j 

X (4(k,7/)4*(k,7/0>, 

where jo(x) = sinx/x is a spherical Bessel function. 

We can factor out the temporal dependence of $ for each k-mode, 

IKM) = £ (k)/kfa) With /k(^r) = 0, /k(7/o) = l- (2-10) 

We will assume that the present day (77 = rj ) gravitational potential field is statistically homogeneous 
and isotropic and hence the expectation of the product, $ (k)$o(k'), is given by the gravitational 
potential power spectrum today, which we call Q; since Q can only depend on k = |k|, we have 

($ (k)$*(k')> = (2nfQ(k) *( 3 >(k - k') » VQ(k)5 kk , . (2.11) 

Thus, the correlation function may be written, 

f d 3 k f Vo f Vo I 

C LT (a) = 4 J j^-^Q(k) J drj J drj' jo(k^ (rj - r]) 2 + (rj - r]') 2 - 2(r] - r])(rj - 77') cosa) 

x h(v)'Mv') ■ 

(2-12) 

Since the mean square anisotropy is just the sum of the anisotropies from the different k modes, 
we may optimize the different k modes (to give minimal C(a)) independently. Also note that the 
optimization just depends on k = |k| and not on the direction of k, so the solution for the /k with the 



same k will be exactly the same. Hence we can assume that only depends on k and we henceforth 
use the notation f k . 

By contrast, for the usual Sachs- Wolfe effect (i.e., with primordial adiabatic fluctuations, and 
purely gravitational evolution thereafter), the temperature correlation function is instead 

1 f cfik 

Csw(a) = g J Q(k)j (2k( Vo - rj r ) sin(a/2)) . (2.13) 

To reiterate: for primordial, linear adiabatic fluctuations evolving purely gravitationally in an O = 
1, matter-dominated universe, /k(^) = and (2.12) vanishes; in this case, (2.13) gives the entire 
anisotropy on large angular scales. Here, we are considering the 'opposite' case in which the primordial 
anisotropy Csw( a ) = 0, and we are seeking to minimize the anisotropy (2.12) arising from the time- 
dependent gravitational potential, independent of any assumptions about gravitational evolution. 

To make contact with observations of large-scale structure (see §IV), it is useful to relate Q(k) 
to P(k), the power spectrum of density fluctuations. Defining the Fourier transform 5 k of the density 
field 5p(x, t])/p as in (2.5), the density power spectrum is defined by analogy with (2.11), 

(5 k 5l,) = (2n) 3 P(k)5 3 (k-k>). (2.14) 

By Poisson's equation, V 2 $ = AixGp, in an O = 1 universe, 

Q(k) = l§ p (k), J = 9 Sj dk k ~ 2p(k) - (2 - 15) 



Minimizing C(0) 

We first consider the minimization of C(0) = ((AT/T) 2 ), the rms variance in the temperature 
fluctuation on the sky; although this is an unmeasurable quantity, it generally sets the scale for the 
temperature perturbations for a model. In the late-time scenario, from (2.12) and (2.15), 

f d 3 k P n° 

C LT (0)=4j ^Q(k) J dr, J dr,'j (k(r,'- V ))f k ( V )f k (r,') 



2tt 2 
_9H$ 



fdkk- 2 P(k) [ n °d V [° drfjo(ktf-ri))f k (v)htf) (2-16) 

J dkk- 2 P(k)I k [f k ] , 



~ 2tt 2 

while, for primordial adiabatic perturbations, the Sachs- Wolfe result (2.13) is 

Csw(0) = - / j^Q(k) = -4 / dkk- 2 P(k). (2.17) 



9 J (2tt) 3 ^ w 8tt 2 

In general, for power spectra that behave as P(k) oc k n for small k, the integrals (2.16) and (2.17) 
will diverge at long wavelengths if n < 2. Since C(0) is not an observable, this divergence is not 
problematic. For example, we would obtain finite results if we calculate physical quantities such 



as a/2(C(0) — C(a)), the rms temperature difference measured by a two-beam experiment with a 
beam-throw of angle a. (In addition, we should take into account the finite width of any real beam.) 
Alternately, we could subtract off the unmeasurable monopole and dipole terms from C(a) (see 
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below) . Here, we are interested in comparing the usual Sachs- Wolfe and the late-time anisotropics in 
/c-space, so we need not perform the divergent wavenumber integrals. 

If we wish to minimize the functional (2.16), it remains to find a set of optimal solutions fk(v)i 
given by the minimization of 

h[fk}= J drj J drfgk(v-rf)'fk(v)fkW)i where 9k{v~v) = ~^~TZ ^ ' 9k(%) = 9k(~x) , 

(2.18) 

with respect to the function fk, which must satisfy the boundary conditions (2.10). 
High and Low Frequency Limits: The Linear Model 

To gain some insight into the level of minimal anisotropy expected, we can analytically explore 
the integral Ik in the short and long wavelength limits, which correspond roughly to krj much larger 
and smaller than one respectively. (It is useful to recall the conversion k = (6.67 x 10- 5 krj o ) h Mpc" 1 .) 
If the the function / is well approximated by its Taylor series close to every point, then in the short 
wavelength limit k — > oo, we may make the substitution in the integral (2.18) 

9k(v - rf) = ™ k ff_~f - *6(k( V - rf)) . (2.19) 

In this limit, Ik reduces to the action functional for a free particle with 'coordinate' fk, and the 
condition for an extremal time history is 

/ fc = 0, i.e., fk(v) = fk n =^HL^IL. (2.20) 

Vo - Vr Vo 

In this case, the minimal temperature correlation function becomes 

c -<°>W(Slr " (2 - 21) 

where Vd = Vo — Vr ~ Vo- Comparing with the usual Sachs- Wolfe expression (2.17), we see that, for 
the same amplitude of present structure P(k), the minimal late-time anisotropy is only smaller than 
the primordial anisotropy in the wavenumber range 

ki] d > 36tt or A < — . (2.22a) 

18 

In the spatially flat cosmology we have assumed, this corresponds to 

AT 1 < 53/i _1 Mpc or A < 333/i _1 Mpc. (2.22b) 

corresponding to a few degrees on the sky. (This lengthscale is comparable to the largest scales 
currently probed by redshift surveys, Cf. Fig. 3 below.) Thus, for the very large-scale perturbations 
which are now starting to be probed by COBE and other experiments, one cannot really do better 
than primordial adiabatic perturbations in minimizing the anisotropy. We can understand results 
(2.21-22) heuristically, by considering the contribution from perturbations of comoving wavelength 
A to the rms anisotropy. For primordial adiabatic perturbations, (AT/T)$w = (1/3)$, where the 
potential fluctuation on scale A is $^ ~ (Sp/ p)\(X/to) 2 - For late-time perturbations, on the other 
hand, the anisotropy is proportional to the integrated time-derivative of the potential, (AT/T)lt ~ 



7 



f&\dt ~ $o,a(V*o)-^^ 2 ) where A^ ~ to/X is roughly the number of lumps of size A between the 
observer and the hyper surface when the potential began to increase (assumed to be at z > 1 here). 
As a result, we find (AT/T)^ ~ 3 ) o,a(-V^o) 1 ' / ' 2 ; n °te that this wavelength dependence agrees with 
(2.21). Comparing this late-time expression with the Sachs- Wolfe anisotropy above, we see that the 
minimal late-time anisotropy is smaller than the primordial anisotropy only for small wavelengths, 
(X/Hq 1 ) 1 / 2 < 1/3, in agreement with (2.22). At large wavelengths, the minimal late-time anisotropy 
is greater than the standard Sachs- Wolfe result because, for primordial adiabatic perturbations, the 
anisotropy is only 1/3 of the gravitational potential fluctuation. This famous 1/3 factor arises from 
a partial cancellation between the gravitational redshift, $, and the varying radiation temperature 
at recombination (AT/T) 7 = —(2/3)$ (for superhorizon perturbations in Newtonian gauge). In 
contrast, for the late-time scenario there is no corresponding cancellation since the last-scattering 
surface is unperturbed. 

The linear form (2.20), f^ n (r]) = (v~Vr)/Vd: is a good paradigm for a slowly increasing potential 
fluctuation, and we shall use it as a fiducial reference against which to compare other results. Let us 
see what it gives for arbitrary k. In this case we have 

^ = J! * J! ^EH* ~f 

Jo Jo K VdV x - U) ^ 23) 

kr] d Si(ki] d ) - (1 - cos(kr] d )y 

(kVd) 2 

where Si is the sine integral function (see Appendix A). Since Si(oo) = ^ and Si(x) = x + 0(x 3 ) for 
small x, we find the limits 

4 in -> ^- for krj d -> oo and /| n -> 1 for krj d -> 0. (2.24) 

The ki] d — > oo limit here agrees with eqn.(2.21). 
Angular Structure for the Linear Model 

Although the absolute values of the temperature fluctuations in late-time models may be com- 
parable to those in standard primordial adiabatic perturbation scenarios, the angular dependence of 
the anisotropy for the two cases can be quite different. To see this, we compare the quantities J(/c, a) 
that are integrated with the power spectrum in the temperature correlation function: 

r4 



C(0) - C(a) = ^ J dk k- 2 P(k)J(k, a) 
J sw (k,a) = - [1 - j (2/c?7osin(a/2))] 



J L t(&, a) =4 dr] dr\ 

Jr/r Jri r 



jo(k(v - v')) -jo(k\ (vo - V) 2 + (Vo - V') 2 - 2(Vo ~ v)(Vo ~ V 1 ) cos a) 



(2.25) 

where we have used the linear model f^ n {rf) in Jlt- In Fig. 1, we plot Jlt and Jsw as functions of 
k for different values of the angle a. Two differences between them are immediately apparent. First, 
consider the behavior at large k: in this regime, the contribution to J is dominated by the a = 
part, J(k,a) — > For the Sachs- Wolfe case, with = 1/36, the contribution does not fall off at 
these small scales. For the late-time model, J(/c, a) ~ 41^ — > 4n/kr] for large k (for a particular value 
of the angle a, this limit is appropriate for all k beyond the maximum of Jlt for that angle — i.e., 
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beyond the scale that contributes the greatest to anisotropies of that angular separation.) Because 
of this difference, power spectra with significant small-scale power may imply greater anisotropies for 
primordial perturbations than in the late-time case. 

Of greater import, however, are the values of J on intermediate scales. Although the maxima of 
J are located at similar values of k for both cases, the maximal values of J are much greater for the 
late-time scenario. On very small angles, the contribution of intermediate scales to the primordial 
Sachs-Wolfe effect is larger than for the late-time scenario; however, as the angle increases, the 
contribution to the anisotropy for the late-time model continues to grow, whereas in the Sachs- Wolfe 
case the maximum value of J remains approximately constant, independent of a. For a given power 
spectrum, then, the anisotropy at a given angular scale will generically be larger for a late-time 
scenario than for primordial perturbations, unless the evolution of the potential is specifically chosen 
to decrease the anisotropy on that scale. 

Low Frequencies 

In the other extreme of low frequencies, or long wavelengths, from (2.18) we have g\. — > 1, so the 
dependence of 1^ on the functions fk drops out; from the boundary conditions (2.10), we then find 
Ik — > 1, which agrees with and generalizes the low frequency result for I k m in (2.24). In this limit, 
from (2.16) we have 



which is 36 times the value for primordial adiabatic perturbations, Csw(0), in eqn. (2.17). That 
is, in the long- wavelength limit, we expect the rms anisotropy to be 6 times larger for the late-time 
scenario than for primordial adiabatic fluctuations; it is worth noting that, for primordial isocurvature 
fluctuations, the rms Sachs- Wolfe anisotropy is also six times larger than for adiabatic perturbations, 



Based on this discussion, one may be tempted to conclude that the behavior min[ifc] ^ 1 as 
k — > is generic. However this is not the case. To see this, consider minimizing 1^ for the small space 
of functions 




as k — > 0, 



(2.26) 



(AT/T) sw>isoc = 2$. 




(2.27) 



In this case, Ik is given by 



h = hn,k + a k h + a k h. 



(2.28) 



where I\ m ^ is given by (2.23), and (see Appendix A) 



h =2tt 




^Ci(|^ d + 27r|)-Ci(|^ d -27r|) 
-^-(kVd) 2 + 0((k Vd ) 4 ) as 



(cos(27rir) + cos(27ry)) 



In 



krj d — 2n 



kr] d + 2-7T 



) 



(2.29) 



k 







9 



and 

I 2 =(27r) 2 t dx f 1 dy sin ( kr ld(x-y)) CQs{2nx) cos{2 ny) 
Jo Jo kr] d {x-y) 



n 



km 



Ci(|fc?7 d + 2tt|) - Ci(|A;?7 d - 2tt|) - In 



2tt 2 (27i) 2 

+ hf d (Si(/c??d + 2yr) + SKkVd ~ 2yr)) + (^q^p (1 " cos (^)) 

(/o? d ) 4 + 0((£;?7 d ) 6 ) as fc -> 



(2.30) 



20tt 2 

Minimizing with respect to gives 



2/ 2 3(fc77 d ) 

and hence 



a k = -A _> ^ Q7r + ^((fer/d) ) as fc -> 0, (2.31) 



uim[4] = J Un , fc - A- _i + (9((A;7/ d ) 2 ) as fc^O. (2.32) 

a 4i2 y 



The small- fc limit of the minimizing integral Ik is 4/9 of the value estimated above, so the corre- 
sponding minimal C(0) is 4/9 of the value given in (2.26). Thus, for the class of functions (2.27), 
in the long wavelength limit the minimal late-time rms anisotropy is 4 rather than 6 times larger 
than for primordial adiabatic fluctuations. However, in order to achieve this limit, the coefficient 
diverges, and therefore the function fa becomes unbounded, as k — > 0. This pathological behavior is 
not what one would expect for the gravitational potential evolution in a 'realistic' late-time scenario, 
but it provides a lower bound on the anisotropy in the long wavelength limit. One can generalize 
this procedure by adding m terms of the form , m sin(7rm(?7 — r} r )/r}d) to fa in eqn. (2.27); we shall 
use this technique below in our numerical work. Extending the sum to larger m further reduces the 
small-/c limit of Ik from the value we found for m = 2. Also note that in the short wavelength limit, 
k — > 00, cik ; 2 —> and we retrieve the linear solution (2.20). However, as m is increased, the ak ; m fall 
off more slowly with increasing k. 

The lesson we draw from this example is that the minimal late-time anisotropy can be sub- 
stantially smaller than that in the linear model f k m which we have been focusing on, but that this 
reduction is achieved at the cost of introducing a potential function fk(v) which varies rather wildly 
with conformal time. 

Multipole Expansion 

It is often convenient to expand the temperature correlation function in angular multipoles, 

,22. op _i_ 1 r 1 
C{a) = V— —Pe(cosa)C £ where C t = 2tt / dcos 9 P^cos 0)C(0), (2.33) 

— ' 47T / 1 

1=0 J 1 

and the P^'s are Legendre polynomials. For the late-time scenario, the angular integral in (2.33) 
decouples the two integrals over conformal time in Eq.(2.12), and the resulting angular power spectrum 
is 

rVo 



Ce,LT = ~~~ I dkk- 2 P(k) ( ° d V j e (k( Vo - v ))f k (v) 



(2.34) 
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where ji(x) is the spherical Bessel function. Again, we can compare to the Sachs- Wolfe anisotropy 
for primordial, adiabatic perturbations, 

Q,sw = § / dkk- 2 P(k)\j e (krh)] 2 . (2.35) 

The angular spectrum (2.34) points to two curious features of the minimization of Clt(O) which 
was missed above. First, it is apparent that by judicious choice of the functions fk(v) we can make any 
particular multipole moment Cg vanish. (We will see an example of this in §IV below.) Once a given 
multipole is set to zero, this specifies all the and the other multipoles will in general be non-zero. 
However, consider the unphysical case in which the potential turns on instantaneously at conformal 
time 77 = rji, i.e., fk(v) = 0(v ~ Vi)- Then the bracketed expression in (2.34) becomes j(.(k(r] — rjf)). 
In particular, in the limit r\i — > rj , this expression vanishes for all I > 0. Thus, if the potential turned 
on instantaneously at the present time, we have Clt(0) = 36Csw(0), in agreement with (2.26), but 
the anisotropy is hidden in the unobservable monopole term Co. (This possibility did not appear 
in the high and low frequency limit discussion of C(0) above, because we did not subtract off the 
monopole term there.) While a useful theoretical foil, this example is not of direct physical interest: 
an instantaneous turn-on of the potential violates causality. More generally, structure on a given 
scale cannot be made in less than a light-crossing time for that scale. Furthermore, we know that 
structure existed before the present time: conservatively, the gravitational potential corresponding to 
the growing mode density fluctuation was essentially in place by a redshift Z{ > 3, corresponding to 
rjf < O.5?7o- This constraint, which we will impose in calculating observables below, implies that the 
i 7^ multipoles will be non-zero, although the higher multipoles may be relatively suppressed. 

To see this, consider the instantaneous turn-on of the potential at Zf. Then the contribution of 
the fc-mode waveband to the £th moment, relative to that for the primordial Sachs- Wolfe anisotropy, 
is 



dCe jLT /dlnk _ 36 
dCifiw /dink 



jliHVo - m)) 



36 



■if 



l-(l + z f )- 1/2 as k^O, (2.36) 



for I 7^ 0. For example, with Zf = 3, in the long wavelength limit, the late-time quadrupole is larger 
than that for primordial perturbations, but the octopole and higher moments are smaller. On the 
other hand, if the potential turns on rapidly at Z{ >> 1, as would be expected in most plausible 
late-time models, then C^lt = 36C^sw 5 independent of k: in this case, the multipole structure in 
the late-time and primordial scenarios is identical, but the late-time anisotropy is 6 times larger. 

III. Minimization for Finite Beams: Formalism and Numerical Results 

In the previous section, we studied the temperature correlation function under the assumption 
of an infinitesimally small beamwidth. We now wish to consider the expectation of the mean square 
anisotropy of a realistic beam configuration when averaged over the sky and averaged over all ob- 
servers. This is given by some rotationally invariant quadratic moment of the temperature field, 
which, using (2.4), may be written 

f d 2 m /AT AT \ If 1 

T = I — — W(cosa) ( -— (n) -— (m) ) = - / d cos a W (cos a) C (a) 
J 4yr \ T T / Xo ft 2J_ X 

=a S wf Q{k) C dv C dv> /.J w{y) Mk ^ {v ' ~ 770)2 +{v ~ 110)2 - 2y{v '~ vo) {v - vo)) 

x h(v)fk(v')- 

(3.1) 
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Here, as before, a is the angle between n and m, and Mucosa) is a weighting function which is 
determined by the beam configuration. For the temperature correlation function C(a) we have used 
(2.12), except that the upper limits on the rj and rj integrals are r\{ instead of rj . This replacement 
incorporates the assumption that the potential perturbations were constant from the epoch rjf until 
the present, rj , that is, that <£>(??) = for r] > rjf. This corresponds to the statement that the growing 
mode density fluctuations were in place at some minimum redshift Z{. 

Our task is to find the functions fk(v) which minimize T. Following eqns. (2.16, 2.18), we use 
the somewhat more compact notation: 



f r/ 3 k f Vf f Vf 

T=4 Q(k)I k h[fk] = dr, dv'gki^^hWhW) 

9k(v, r/) =^ J _dy W(y) jo (k\J (rf - r) Q ) 2 + (77 - t] q ) 2 -2y(rf- r] Q ) (r] - rj^j 

Note that from the definition of g k , we have g k (r), rj') = g k (rj', r i)i an d in any case the integral I k only 
depends on the symmetric part of g k . 

Explicitly, the condition that f k gives an extremum of the functional I k is that 

^-Iklfk + eA] =0 (3.3) 

u€ e=0 

for all functions A(r]) which are zero at the endpoints, A(r] r ) = A(??f) = 0. Thus 



d 
de 



rm rm 

I k [f k +eA] =2/ drj d V 'g k ( V ,rj')f k ( V )A(rj') 
e=o J Vr J Vi 

rm rm . rm 

= -2/ dr] drj' g kj2 (v,v')fk(v)^(v') + 2 / drjf k (rj)[A(rj { )g k (rj,rj { )-A(rj r )g k (rj,rj T )} 

rm rm 

=2 drj d V 'g k ,2(r]'-v)fk(v)Mv') 

(3.4) 

where we use the notation g k 2 to mean differentiation with respect to the second argument, g kj2 = 
dg(x, y)/dy, and we have used the symmetry of g k . Since this is true for all variations A, we see that 
this is equivalent to the condition 

rm . rm 

/ dr]g kj2 (ViV') fk(v) = 9k,2(VhV') ~ drjg k)l2 {rj, rf) f k (v) = V rj' E fa, r] { ). (3.5) 



This expression is of the form {operator) x f k =function. If we can find the inverse of this linear integral 
operator, then we can solve for / mi \, which minimizes I k . Using this equation in the definition (3.2) 
of I k , integrating by parts, and using the boundary conditions (2.10), we find that for a true extremal 
function 

rm 

min[4] = g k (m,Vf) ~ / drj g k ,i(v,Vf) f mm k(v)- (3-6) 

The details of the numerical procedure we use to find f mm k are given in Appendix B. We replace 
the continuous conformal time interval rj G (i] r , rj{) by a grid of N points rj[, and use the trapezoidal 
rule approximation to convert the double integral for I k into a double sum. Following the discussion 
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of §11, we choose rjf = 0.5r) o and to excellent approximation set r] r = 0. We generally find that the set 
of / mi \(?7i) is discontinuous at the scale of the grid. To impose a smoothness cutoff independent of 
the grid size, we therefore expand the fk in sine waves, generalizing (2.27), 



fk(v) = — + Y] ak,m sin ( mix— ] , (3.7) 

with integer m. As required, this satisfies the boundary conditions fk(0) = 0, fk(vf) = 1- By taking 
n < N, we restrict fk(v) from varying significantly over the grid scale. For the results shown below, 
we used n = 70 or n = 100 sine waves and N = 200 grid points. As a check, we have also expanded the 
fk(v) as polynomials of order n (with appropriate boundary conditions) and minimized the Ik with 
respect to the coefficients using an n — 1-dimensional simplex algorithm. The results thus obtained are 
consistent with the analytic and numerical techniques discussed above (although this "brute force" 
method generally finds Ik larger than the more rigorous discretization). 

Before we can proceed to minimize Ik, we must choose the particular weighting, W, for which 
we wish to minimize T. The functions fk which minimize T are unlikely to be the same for different 
functions W(cosa). First, to make contact with the discussion of §11, consider minimizing the mean 
square anisotropy at a point, C(0) = ((AT/T) 2 ). In this case, the appropriate window function is 
W(y) = lim e ^ + 2S(y — 1 + e), where the small positive e guarantees that the full 5-function is included 
in the integral. Substituting this into (3.2), we find gkiv^v') = sm(k(rj — r]'))/k(rj — 77'), in agreement 
with (2.18). 

To minimize other quantities more directly related to experimental observations, it is convenient 
to again expand in multipoles. Using the multipole expansion (2.33), we can write 

T = J2^r WeCe (3 - 8) 

1=0 



where 



00 1 /"^ 

W(x) = V(2£ + 1) W e P t {x) and W t = - I W(x) P e (x) dx. (3.9) 
T^L 2 J-i 



In terms of we may rewrite the expression for gk in (3.2): 



1 00 «i 
9kM)=-Y^ + ^)Wi\ dyP e (y) 



sin/c A /(?7 / - r] ) 2 + (rj - r] ) 2 - 2 y {rf - r) Q ) (rj - rfaj 



00 



kyj (rf - T]o) 2 + Jji - T] ) 2 -2y(rf- r) ) (rj - rj Q ) 



=\ E( 2 ^ + !) W t I 77 ,_ n v _ =T J t+lW ~ J t+^KV ~ Vo)) (3.10) 
e=o ^Vw VojyV Vo) 

00 

= J2( 2i + 1) WdtW - vo)) MHv - vo)) 

£=0 



where the J n are Bessel functions, and ji(x) = a/tt/ (2x) J i+ i(x) are spherical Bessel functions. 

Minimization for a finite beam experiment 

Although we discussed the minimization of the rms temperature anisotropy C(0) in §11, we also 
pointed out that it is not an observable quantity. Instead, one measures the temperature fluctuation 



13 



over some finite region of the sky determined by the beam pattern of the instrument. In many cases, 
the instrument beam is roughly of Gaussian form. Given the intrinsic temperature pattern on the 
sky, one can construct the sky temperature pattern convolved with the beam. For a Gaussian beam 
of width a (radians), the two-point correlation function of this beam-smoothed pattern is given by 



in (3.10) to determine gkiVirf)- 

As an example which we will use below, consider the recent COBE DMR observations (Smoot, 
etal. 1992). The DMR beam is approximately Gaussian with a 7° FWHM; since a = 0.43 FWHM, 
this implies &COBE = 5.2 x 10~ 2 . The DMR team published three results of interest: the quadrupole 
anisotropy; the correlation function (3.11) with terms i = 0, 1, 2 removed; and the rms fluctuations 
smoothed on 10°, with the monopole and dipole removed. The latter result is the most useful for 
us, and we can express it via (3.11) as C(0, a[10°]), where the beam-width corresponding to the 10° 
FWHM is a [10°] = y/2cr C OBE, and we remove the I = 0, 1 terms from (3.11). 

In §IV, we will give results for C(0, a) for a variety of beamwidths a; from (2.34), (3.11), and 
(3.12), we have 



In Fig. 2, we show the minimizing integral min[ifc(0, a)] for C(0, c[l ]) and C(0, cx[10 ]) as a function 
of krjo, for the late-time scenario (minimized according to eqn. 3.7, and shown as the points denoted 
LT in the figure) and for primordial adiabatic perturbations (curves denoted SW). This shows that, for 
a given power spectrum P(k), the minimal late-time anisotropy on 10° is smaller than for primordial 
adiabatic perturbations, unless the spectrum were narrowly peaked around krjo ~ 20, in which case 
they would be similar. On 1°, the situation from Fig. 2 is less obvious: if P(k) has little power at 
krjo < 20, then the minimal late-time anisotropy on this scale could be larger than the primordial 
anisotropy, but if there is significant power on these large scales, the larger Sachs- Wolfe 1^ in this 
small-/c region would lead to larger relative anisotropy for primordial fluctuations. Note that the 
corresponding potential functions fk(j]) for the minimizing late-time scenarios are wildly oscillating 
functions of r\ (see Appendix B and Fig. 6). If we instead constrained fkiv) to be a more gently 
varying function as in §11 (e.g., eqn. 2.20), the resulting I^s in the late-time models would lie above 
the Sachs- Wolfe results over the range of k shown in Fig. 2 (Cf. Fig. 1). 

Minimization for a Multiple Beam Switching Experiment 

A common type of MBR anisotropy experiment is a switching experiment which, in its simplest 
form, consists of measuring the temperature convolved with a Gaussian beam at 2 or 3 evenly spaced 
points on the sky. For reference, we give here the corresponding window functions. For a 2-beam 




(3.11) 



Thus, to minimize T = C(a, a) we would use 



W e = e-^ £+1 ^ 2 P e (cosa) 



(3.12) 




(3.13) 
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experiment with Gaussian beam-width a, and beam throw a, one should take 
and hence 

W e = e- £ ^ G \l - P e (cosa)). (3.15) 
For a three-beam experiment with throw a between adjacent beams, one can take 

1 / /AT AT AT \ 2 \ 4 1 

r = - ( - 2— 2 + — J j = C(0, a) - -C(a, a) + -C(2a, a), (3.16) 

and hence 

W t = e~ e ^\l - ^P £ (cosa) + ±P e (cos2a)). (3.17) 

The number quoted for "AT/T" for these experiments is y/T times some factor, where the factor 
used may vary from experiment to experiment. 

IV. Power Spectra and Results 

At present, the COBE DMR results are the only probe of the power spectrum on scales larger 
than a few hundred h~ 1 Mpc. Under the standard hypothesis of primordial adiabatic perturbations, 
COBE provides direct information on the large-scale primordial power spectrum through the Sachs- 
Wolfe effect, e.g., eqn.(2.13), and Smoot, etal. (1992) find P(k) oc k n , with n ~ 1 ± 0.5. However, if 
we discard the assumption of primordial perturbations, the results of a MBR anisotropy experiment 
can no longer be used to determine P(k) in the absence of a specific model for the evolution of 
the gravitational potential. Specifically, the consistency of the COBE results with the inflationary 
prediction of a Harrison-Zel'dovich spectrum, P(k) oc k on large scales, could be an artifact of some 
other power spectrum along with suitably chosen /c-dependence for the evolution functions fkiji) (Cf. 
Fig. 1). 

Complementing the COBE results, there have recently been several determinations of the galaxy 
power spectrum from catalogs derived from the IRAS survey and others. However, these observations 
do not determine the power spectrum on scales large enough to overlap those probed by COBE. In 
particular, while COBE probes the (primordial) shape of the spectrum on very large scales, galaxy 
observations only extend up to scales of order 100 h _1 Mpc, where significant processing of the 
primordial spectrum has taken place. If we do interpret the COBE results as a Sachs- Wolfe probe 
of primordial perturbations, then the resulting COBE spectrum (e.g., Harrison-Zel'dovich) need only 
be matched onto the smaller scale galaxy observations (modulo such crucial factors as biasing and 
selection effects). In the late-time scenario, however, the shape of the power spectrum on large scales 
is not uniquely fixed, but only constrained, by the COBE observations. 

In the previous section, we compared the contributions to the anisotropies for the late-time and 
primordial scenarios from a given wavenumber k, for the same amount of power P(k). Here, we 
integrate these contributions over several phenomenological power spectra to make predictions for 
the observable anisotropy for different beam configurations. These phenomenological spectra include 
two models which approach the Harrison-Zel'dovich form at large scales, but which differ on small 
scales: one is the standard cold dark matter (CDM) spectrum, and the other is an analytic fit to the 
QDOT galaxy power spectrum. However, since we are not assuming primordial perturbations, we 
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should not interpret the COBE results to mean that the spectrum must approach something like the 
Harrison-Zel'dovich form at large scales. Therefore, we also consider a third spectrum, based on a 
fit to the QDOT data at small scales as well, but which is more sharply cut-off at large scales, with 
P(k) oc k A as k -> 0. 

The three spectra we use are 

( kT 2 (k) CDM 
P(k) oc I ke- k P/ (1 + (k/k ) 2 ) MGSS-I (4.1) 
[ k 4 e~ k P/ (1 + (k/k ) 5 ) MGSS-II 

Here T(k) is the CDM transfer function of Bond and Efstathiou (1984) for Q = 1, Q B = 0.03, and 
h = 0.5, 

T(k)= 1 + (ak + (bkf /2 + (ck) 2 Y lf \ (4.2) 

where 

a = 5.8(ft/i 2 ) -1 Mpc, b = 2.9(0/i 2 )" 1 Mpc, c = 1.6(0/i 2 )" 1 Mpc, u = 1.25, (4.3) 

and p = 8/i _1 Mpc, k^ 1 = 30/i -1 Mpc. The last two spectra in this list were used by Martinez- 
Gonzalez, Sanz, and Silk (1992) as approximate phenomenological fits to the power spectrum from 
the QDOT survey of IRAS galaxies (Kaiser et al. 1991). These models, along with the power 
spectra inferred from the QDOT (Feldman, Kaiser, and Peacock, in preparation) and 1.2 Jansky 
IRAS redshift catalogs (Fisher et al. 1992), are shown in Fig. 3. We normalize the model spectra in 
the usual way, by setting the rms mass fluctuation within spheres of radius 8 h _1 Mpc to be 

4 = li^f) ) = ^2 [°° dkk 2 P(k)W 2 (kR)\ R=8h -i Mpc , (4.4) 

where the window function is W(kR) = 3(sin/cR — kRcos kR)/(kR) s . Below, we present results for 
Ta 8 2 . 

In Fig. 4, we show the anisotropy expected for a 2-beam experiment at a given angular scale, 
(AT/T) 2 , /2 = (C(0) - C(0)) in units of of, for the linear late-time model f^ n of eqn.(2.20) and 
for primordial adiabatic perturbations, for the three spectra of eqn. (4.1). For the MGSS-1 and 
CDM spectra, which both approach the Harrison-Zel'dovich form at large scales, the large-angle rms 
anisotropy (AT/T) e in the late-time model is roughly three times larger than the corresponding 
anisotropy in the model with primordial adiabatic perturbations. The results are different for the 
sharply falling spectrum MGSS-II: since this spectrum has no power on large scales, the late-time 
anisotropy is smaller than the primordial Sachs- Wolfe anisotropy in this case, but by less than a 
factor of 2; this is in agreement with the expectation from eqn. (2.22). Note also that, in accord with 
Fig. 1, for the same power spectrum P(k), the angular dependence of the anisotropy for the linear 
late-time model differs substantially from that for primordial fluctuations at small angles. 

In Figure 5, we show the correlation function C(0, a) as a function of beam- width a, with the 
monopole and dipole terms subtracted off. Again we use the spectra of (4.1), and show results for 
primordial adiabatic perturbations (eqns. 2.35, 3.11, and 3.12) and for late-time perturbations, eqn. 
(3.13), minimized according to eqn.(3.7). A note of caution in reading the late-time curves in this 
Figure: the integral ifc(0, a) has been minimized independently at each value of a, i.e., different 
potential functions fk(v) have been chosen at each a. Therefore, a given late-time curve in this figure 
does not correspond to a fixed late-time scenario (i.e., to a fixed set of fk(v))i but rather to many 
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different scenarios. Consequently, the cr-dependence of the late-time curves should not be interpreted 
as implying that the anisotropy for a given late-time scenario falls off with increasing beam-width 
according to these curves. In fact, for a fixed late-time model that minimizes the signal at some 
particular <r c , the fall-off at a > a c would be more gradual than in the figure, while the rise in the 
signal at a < a c would be steeper. 

The COBE DMR result for the fluctuation on 10 degrees, C DM r(0, a[10°]) = (1.2 ±0.4) x 10~ 10 , 
is shown for comparison (Smoot et al. 1992). For beam- widths less than a few degrees, the minimal 
late-time anisotropy is comparable to the primordial adiabatic result. However, at larger beam-widths, 
the minimal late-time result falls sharply below the Sachs- Wolfe anisotropy. We can understand this 
result heuristically as follows. For very large beamwidth a, the contribution of higher multipoles to 
the sum in (3.13) is strongly suppressed. As a result, in this limit, the anisotropy is dominated by 
the quadrupole (and perhaps the octopole). However, from the discussion following (2.34), it is clear 
that one can choose a set of functions fk(v) to make a particular multipole Cg, e.g., the quadrupole, 
vanish. For r/f sufficiently small, one would normally expect this to produce large values for the other 
multipoles, but these higher moments enter the large-cx anisotropy with very small weighting. As a 
consequence, the minimal anisotropy for large beam-widths can be quite small. 
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V. Conclusion 

We have seen that the MBR anisotropy signature of post-recombination structure formation gen- 
erally differs from that of primordial fluctuations. Consequently, once MBR anisotropy experiments 
and large-scale structure observations begin to overlap significantly in the lengthscales they probe, 
comparison of the two would allow one to definitively test whether the fluctuations are primordial or 
more recent in origin. We have found that, for a given amplitude of present large-scale structure P(k), 
the minimal late-time anisotropy can be up to an order of magnitude smaller than the correspond- 
ing Sachs- Wolfe anisotropy for primordial adiabatic perturbations. However, as comparison of Figs. 
4,5, and 6 show, this minimum is only achieved if we allow sufficiently pathological time dependence 
for the gravitational potential $>k(v)- If* we restrict the time-dependence of the potential to more 
well-behaved forms more plausibly to be expected in late-time scenarios (e.g., the linear model or the 
high-redshift, rapid turn-on model of §11), then the large-angle anisotropy in the late-time scenario is 
generally comparable to or larger than that due to primordial fluctuations. In particular, this will be 
the case if the present density fluctuations have substantial power on scales larger than k' 1 ~ 53/i _1 
Mpc (Cf. eqn. 2.22). While this result runs counter to part of the motivation for late-time phase 
transitions, it is not necessarily a negative result for them, given that large-scale MBR anisotropies 
have now been observed. Furthermore, as recent work suggests (Frieman, Watkins, and Hill, in 
preparation), perhaps the most likely role for late-time transitions is to amplify perturbations that 
were initially present over some range of wavelength, rather than to replace primordial fluctuations 
entirely. In this case, the final power spectrum is due to a combination of primordial and late-time 
effects, and the induced anisotropy will correspondingly arise from both. 

We should also comment on the relation of our work to the recent paper by Martinez-Gonzalez, 
Sanz, and Silk (1992). These authors calculate the contribution from the time- varying gravitational 
potential to the anisotropy, as do we. However, they considered a very specific mechanism for the 
time evolution, namely, that the potential varies in time due to mild nonlinearity of the density 
inhomogeneities; this effect occurs even in an Einstein-de Sitter (spatially flat) cosmology. The 
anisotropy induced from this non-linear gravitational evolution is small, 5T/T ~ 10~ 6 , and can be 
more than an order of magnitude below that arising from primordial adiabatic fluctuations. This 
small number is not to be compared with ours, since to the second order effect they have calculated 
must be added either: 1) the effects associated with growing the perturbation to the amplitude at 
which second order effects become important, or 2) the primordial anisotropy from last scattering. 
For the standard gravitational instability scenario, effect (2) dominates over the second order effect 
calculated by Martinez-Gonzalez, etal., unless the universe is reionized after recombination. For late- 
time scenarios in which (2) is absent or negligible, we have shown in this paper that effect (1) is not 
necessarily very small and that, in 'realistic' models, it is likely to dominate over the second order 
gravitational contribution to the anisotropy. 

Finally, it is worth noting that our methods could be extended or applied in a number of 
ways. For example, one could use them to estimate the expected anisotropy in a specific late-time 
scenario, in topological defect models of structure formation, and in the loitering universe model 
(Sahni, Feldman, and Stebbins 1992, Feldman and Evrard 1992). In addition, one could consider 
models in which the perturbations induced at late times are non-Gaussian, which one might expect 
to be a natural outcome of late-time transitions. 
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Appendix A 

We evaluate here several of the integrals in the text. First we compute the angular average in 
eqn. (2.8). Consider two unit vectors n and rh separated by angle a, n • m = cos a, and a wavevector 
k such that n • k = \x = cos (3 and ih • k = sin a sin (3 cos ip + cos a cos (3. Then 



(e 



jk-(an— 6m)\ 



J dcos (3 dip exp [ik • (an — 6m)] 



/n-m=cosQ 



/ dcos P dip 

= — J dcos (3 dip exp ik(a cos/3 — 6 sin a sin/3 cos ip — 6 cos a cos (3) 

= ^- J dcos /3 exp [ik(a-bcosa) cos/3] J dipexp(— ikb sin a sin (3 cos ip) 

= -j— I dcos (3 exp [ih (a — b cos a) cos f3] x 2TiJo(bk sin a sin (3) 
An J 

djJL cos [/ik(a — 6 cos a)] Jq (bk sin a a/ 1 — // 2 ^ 



Doing the remaining integral and applying to the expression (2.8) gives 

/ -ik-n(r7o-??)+«k / -m(77o-r? / ) 



(A.l) 



)n-m=cosa = 3o(k\ (Vo - v) 2 + (Vo ~ v') 2 ~ 2(Vo ~ v)(Vo - rf) cos a), (A.2) 



where jo(x) = sinx/x. 

For the integral /[ m in eqn. (2.23), we have 



rliii 



= [ dx [ 
Jo Jo 



1 



2.7-1 j\v x 
1 

dv(l - \v\) 



sin(ki] d (x - y) 

„ * k Vdix -v)) 

sin(kr] d v) 



2-\v\ 



du 

krjd v 

sin(kr] d v) 



u = x + y, v = x — y 



(A.3) 



= 2 



l krj d v 
kyd Si(fc??d) ~ (1 ~ cos(kr] d )) 
(km) 2 
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Next consider the integrals I\ and I2 in eqns.(2.27) and following. We have 

sin(ki] d (x - y)) 



h = 2tt 



[ dx [ 
Jo Jo 



dy 



(cos(27rx) + cos(2ny)) 



2;; / dv 



ki] d (x - y) 

du — ; cos(ttu) cos(nv) u = x + y, v = x — y 



-4 



dv 

1 ki] d v 



krj d v 
cos(7rz;) sin(7r|t;| 



Z 1 , sm(kriAv) . , n 
-4 / dv — , - ' sin(27rw) 



(A.4) 







krj d v 



= 2 



Ci(|/c77 d + 2n\ v) - Ci{\kr] d - 2tt\ v) 
kVd 



v=l 



J v=0 



ki] d 



Ci(\krj d + 2tt|) - Ci(\krj d - 2ir\) - In 



ki] d + 2n 



krj d — 2n 



-^(k Vd ) 2 + 0((k Vd ) 



and 



! f dV f 

J-l J\v 



h = (2tt) 2 dx dy 



= 2tt 



7T 



= 7T 



sin(krj d (x - y)) 



cos(27rx) cos(27ry) 

kr] d (x - y) 

2 ~' u ' sin(/c?7 d t;) 

aw — ; cos(7r(it + V)) cos{2n{u — v)) u = x + y, v = x — y 



1 sin(/c?y d w) 

ai> — : 

1 K?y d w 



/c77 d v 

-- sin(2yr|z;|) + 2(1 - \v\) cos(2ttw) 

TV 



Ci(\krj d + 2tt\v) - Gi(\krj d - 2ir\v) 



Att 



kVd 

2 1 - cos(kr] d ) 



v=l 



+ 2tt 



v=0 



n 



ki] d 



(k Vd ) 2 - (2tt)2 
Ci(|fc + 27r|)-Ci(|fc-27r|)-ln 

(2tt) 2 



k + 2n 



k-2ix 



2 Si(kr] d + 2%) + Si(kt] d — 2tv) 
kVd 



2ir 2 

+ — (Si(Jfe + 2tt) + Si(Jfe - 2tt)) 



(A.5) 



+ 



(2tt) 2 - (^ d )2 



(1 - cos(/c?7 d )) 



20tt 2 



(^ d ) 4 + 0((^ d ) 6 ) 



In the above expressions, we have used the sine and cosine integrals, 

P0 ° cost 



Si(x) = [ X —dt Ci(x) = - f 
Jo t J x 



dt, 



(A.6) 
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Appendix B. Numerical Minimization of the Integral 1^ 

We wish to find the function F(x) = fk(xVf) which minimizes the integral 

1= f dx f dx'G(x,x')F(x)F(x') (B.l) 
Jo Jo 

with the boundary condition that 

F(0) = ^(lHl. (B.2) 

It will be useful to have the symmetry G(x,x') = G(x',x) in this integral. The function G( , ) need 
not be symmetric with respect to its two arguments, but clearly the integral only depends on the 
symmetric part. If G{ , ) is not symmetric, then we can use instead 

G(x,x') -> -(G(x,x') +G(x',x)) (B.3) 
which has the required symmetry. We may integrate (B.l) by parts to obtain 

7=G(1,1)- C dxG,i(x, 1) f(x) - t dx'G, 2 (l,x')f(x')+ f dx t dx 1 G, 12 (x, x') f(x) f(x') 
Jo Jo Jo Jo 

=G(l,l)-2 [ dxG A (x,l) f(x)+ [ dx [ dx'G i 12 (x,x') f(x) f(x') 
Jo Jo Jo 

(B.4) 

where ; i denotes differentiation with respect to the first argument and ^ the second. 
Trapezoidal Rule Discetization 

We can discretize this by only specifying F(x) at a finite ordered set of points {xi} for i = 
0, . . . , N + 1 with xq = and xn + \ = 1. First consider the trapezoidal rule approximation to (B.l): 

N+l N+l 

I=J2J2 °ijFi Ax i F 3 Ax 3 ( B - 5 ) 
i=0 j=0 

where Fi = F(xi), and 

( \xi t = f 2^T(^i-^o) ^ = 

Gij = G(xi,Xj) A Xi = I ±(x i+1 -Xi-i) l<i<N Fi=! 5257(^+1-^-1) 1 < % < N 

( 1(1 -x N ) i = N + l i = N+l 

(B.6) 

Then using the notation 

f Fi - F i = 

AFi = \ (F i+1 - 1 < z < N (B.7) 

[ (F N+1 -F N ) i = N + l 

and the fact that Fo = and F/v+i = 1 to rewrite (B.5) and then "difference by parts" to obtain the 
"difference" analog of (B.4): 



N+l N+l N N N 

1 = 4 E E Gij AFi AFj = YsY. MijFiFj - 2 £ fl^ + C (B.8) 



i=0 j=0 i=l j=l i=l 
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where 

Mij =-(Gi+i j+i - Gi+i j-i - Gi-i j + i + Gi-i j-i) 
Bi =-(G i+ i jv+i + G i+ i n - Gi-i n+1 - Gi-i n) (B-9) 

C =-(Gnn + Gn+i n + Gn n+1 + Gn+i jv+i) 

One may then use standard linear algebra techniques to solve for the set of {Fi} which extremizes 
this discretized approximation to /: 

N 

Fr = J2 M v 1 B J . (B.10) 

i=i 

For large enough N, the set {F™ 11 } should give a good approximation to the function F(x) which 
extremizes I. Note that for {FP 1111 } to actually be a minimum and not a saddle-point requires My to 
have only positive eigenvalues. While it is clear that this is true of the continuous integral operator, 
we have not shown it for this discretized representation. Also note that, instead of the trapezoidal 
rule, we could use some higher order approximation to the integral. As long as this approximation is 
bilinear in Fi, this would just correspond to a different matrix Gij, and the rest of the analysis would 
carry through. 

With the technique described above, we generally find that the set of {F™ 111 } is discontinuous at 
the scale of the grid. The discontinuous nature of the minimal {F- 11111 } holds even for small n, so we 
are confident that this behavior is not a result of round-off error. The convergence of / for increasing 
N is also fairly rapid, although there is significant "noise" in this convergence. We cannot be certain 
that the minimal I we obtain with this method is obtainable with any smooth function F(x). While 
there is a piecewise bilinear 2-dimensional integrand which gives this integral, we cannot be sure it 
is of the form G^2(x, x') F(x) F(x'). We also cannot be certain that My has no negative eigenvalues 
and that the "extremum" we have found is not actually a saddle point. However, we believe that the 
the minimal values we obtain are actual lower limits to what is achievable with a smooth function. 
In any case, most of the peculiarities of these results are almost certainly dependent on the difference 
scheme used. The basic problem is that the integration scheme requires the function to be smooth on 
the grid scale in order to be accurate, while the minimization scheme is forcing the integrand toward 
discontinuity. The solution to this problem is to have two resolution scales, the smaller one used to 
perform the integral and the larger one setting a bound on the jaggedness of F(x). This will allow 
us to vary the accuracy of the integral and the smoothness of F(x) independently. 

Sine Wave Expansion 

One way to implement a smoothness cutoff, independently of the size of the grid, is to first 
expand F(x) in some set of smooth, linearly independent functions. The simplest example are sine 
functions, 

n 

F(x) = x + ^ a a sina7rx (B.ll) 

a=l 

for integer a. Since sina7rx = for x = 0, 1, eq. (B.ll) enforces the boundary condition (B.2). 
The sina7ra; for integral a are a complete set of linearly independent functions on the interval (0, 1), 
and will be linearly independent of the function x for finite n. However, since x and the {sincwx} 
are nearly not linearly independent for large n (i.e., smoothed over an interval of size ~ 1/n, one 
can approximate x by a superposition of n sine waves), one must be careful about unnecessary 
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small eigenvalues in the matrix L ab below, which may cause problems numerically. An expansion in 
Legendre polynomials such as 



F{x) = Pi (X) + J> a (P2a+l(x) - Pl(x)) (B.ll') 
a=l 

would not have this problem. However, sines are superior to Legendre polynomials, because they 
have a fairly uniform variation over the interval and hence give a fairly uniform resolution. The P^'s 
for large / have much more rapid variation near x = 1 than near x = 0. 

Using a grid, {xi}, as above and the sine wave expansion (B.ll), we may rewrite (B.9) as 

n n n 

I = Y.H L abu a u b - 2 Paa a + D (B.12) 

a—l 6=1 a—1 

where 

N N N N N 

Lab = E E M iJ 8 a4 Pa = Yl BiS * "EE Mi 3 XiS a 
i=l j—1 i=l i=l j=l 

N N N 

D = MijXiXj — 2 B{Xi + C s a = sin aixxi 

i=i j=i i=i 

This has extremum 



(B.13) 



„ mm 



= L^ft, F min (x) =x + J2 «f n sinanx. (B.14) 

6=1 a=l 



Of course, we must have n < N for to be non-singular. If n = N, then L a (, is related to M a (, 
by a similarity transformation, and nothing has been changed. However, for n < N we prevent 
F(x) from varying wildly on the grid scale. In Fig. 6, we show the function fk(^]) for the solution 
which minimizes the integral lfc(0, a [10°]) (Cf. 3.13). Here, we have fixed kr/o = 20 and iV = 200 
grid points, and we vary the number of sine waves used, n = 2, 10, 99, and 100. We see that, as 
n is increased, the function becomes increasingly noisy, and that, even at large n, the form of 
the potential function can change substantially with a small increment in the number of sine waves. 
Nevertheless, for sufficiently large n, the value of the integral 1^ converges fairly well. Comparison of 
the integral for large and small n shows that the linear model n = overestimates I\. by more than 
an order of magnitude. 
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Figure Captions 

Fig. 1: The functions J(k,a) (eqn. 2.25) as a function of kr] for the primordial Sachs- Wolfe (SW) 
anisotropy (dotted curves) and for the linear late-time model (LT), f^ n (solid curves). Moving from 
bottom to top, the curves correspond to angles a = 10, 20, 90, and 180 degrees. 

Fig. 2: The integrals Ik for C(0,cr[l°]) and C(0,cr[10°]) for primordial adiabatic perturbations (SW 
curves) and for the late-time scenarios (LT points) which minimize the correlation function filtered 
through these two beam- widths, cf. eqn. (3. 13). Note that the monopole and dipole terms have been 
removed. 

Fig. 3: Density power spectra P(k) (eqn. 4.1) are shown for the three models of eqn.(4.1), denoted 
MGSS-1 (solid curve), MGSS-2 (dot-dash curve), and CDM (dashed curve), all normalized to ag = 1. 
Also shown are the inferred galaxy power spectra from the QDOT (open squares, Kaiser, etal. 1991, 
as reanalyzed by Feldman and Kaiser, in preparation) and 1.2 Jansky (crosses, from Fisher, etal. 
1992) redshift surveys based on the IRAS catalog. The survey spectra have not been corrected for 
redshift distortions, and are shown here principally to motivate the phenomenological fits of eqn. (4.1). 
(Note that our convention results in a factor (-7r/2) 3 difference in the value of P(k) from that used by 
Fisher, etal. (1992) for the 1.2 Jansky results.) 

Fig. 4: The temperature correlation function [C(0) — C(9)]a^ 2 vs. 9 is shown for the linear late-time 
model f l k m of eqn. (2.20) (curves marked by crosses, triangles, and boxes) and for primordial adiabatic 
fluctuations (unadorned curves), for the three phenomenological spectra of eqn. (4.1). 

Fig. 5: The temperature correlation function C(0, a[FWHM])a^ 2 for an experiment of Gaussian 
beamwidth a is plotted as a function of beam FWHM, for primordial adiabatic perturbations (un- 
adorned curves) and for the late-time scenario (curves marked by crosses, triangles, and boxes) nu- 
merically minimized according to eqn. (3.7) (Cf. eqn. 3. 13), again for the three spectra of (4.1). Note 
that each late-time curve corresponds to many different late-time models, each of which minimizes 
the anisotropy at a given a. The COBE observation at FWHM of 10 degrees is shown by the closed 
circle. 

Fig. 6: The potential function fk(v) as a function of conformal time 77/770 is shown for the late-time 
model that minimizes the 10° anisotropy C(0, a [10°]), for fixed wavenumber kr/o = 20, N = 200 grid 
points rji, and for n = 2, 10, 99, and 100 sine waves. This demonstrates that the time dependence of 
the gravitational potential becomes increasingly noisy for large n, but that the integral 1^ converges 
as the number of waves is increased. Note that the boundary conditions fk(0) = 0, fk(f]f = 0.5?7o) = 1 
have been imposed. 
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